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Abstract 

Linear optical functions of cubic and hexagonal BN have been studied within 

first principles DFT-LDA theory. Calculated energy-loss functions compare 

well with experiments and previous theoretical results both for h-BN and for 

c-BN. Discrepancies arise between theoretical results and experiments in the 

imaginary part of the dielectric function for c-BN. Possible explanation to 

this mismatch are proposed and evaluated; lattice constant variations, h-BN 

contamination in c-BN samples and self-energy effects. 
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I. INTRODUCTION 



The properties of boron nitride (BN) have motivated detailed theoretical and experimen- 
tal studies since a long time |jl|J^,^ . Many advanced technologies rely on Boron Nitride and 
on materials based on it, due the wide spectrum of properties offered by its polymorphic 
modifications, two graphite-like and two dense ones. Boron nitride shares many of its proper- 
ties, structures, processings and applications with carbon. Cubic boron nitride (also known 
as sphalerite boron nitride and abbreviated as Z-BN, c-BN or /5-BN),with sp^-hybridized 
B-N bonds, has the diamond crystal structures and a similar lattice constants. Its physical 
properties, such as extreme hardness, wide energy bandgap, low dielectric constant and high 
thermal conductivity, are also very near to those of diamond. These amazing properties of 
c-BN have many appealing applications in modern microeletronic devices, and make it use- 
ful also as a protective coating material or in high-duty tools Hexagonal BN (h-BN or 
a-BN), an sp^-bounded layered compound and graphite also resemble each other in term 
of crystal structures, lattice constants, and physical properties such as strong anisotropy. 
Due to its high termal stability h-BN is a widely used material in vacuum technology. It 
has been employed in microelectronic devices 0, for x-ray lithography masks 0, and as a 
wear-resistant lubrificant [0. The hexagonal phase is also the underlying structure of the 
BN nanotubes, which are systems of growing interest nowadays P,P|JTm. 

On the other hand, significant differences exist between carbon and boron nitride, due 
mainly to differences in their chemical bonding. Infact h-BN and c-BN mechanical strengths, 
thermal conductivities and Debye temperatures are lower than in their carbon counterparts. 
Moreover h-BN is electrically an insulator, while graphite is a conductive semi-metal. In 
contrast to diamond, which can be readly doped only as a p-type, c-BN can be doped either 
n- or p-type. Moreover c-BN does not react with ferrous materials, even at high temperatures 
(T~ 1600/^) [0, and, last but not least, both c- and h-BN are more resistent to oxidation 
than their carbon counterparts 

An important issue is the experimental evidences for negative electron affinity (NEA) in 
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cubic boron nitride samples [0,^. The electron affinity of a semiconductor is the difference 
between the vacuum energy level and the conduction band minimum level. NEA occurs when 
the conduction band minimum lies above the vacuum energy level. Any electron promoted 
into the conduction band has then enough energy to escape into vacuum. One of the most 
striking application is in electron emission devices, to obtain the highest electron emission 
density at the least energy expenditure (cold cathode emitters) 11^. Correct evaluation of 
the one electron transition energies and the position of the single particle states turn out to 
be fundamental points relatively to this issue, and will be here addressed in detail. 

It is now generally agreed experimentally that bulk c-BN is the thermodynamically stable 
phase at ambient conditions, and that the less dense h-BN becomes stable at temperatures 



exceding T ~ 1200/^ [|T5|. Due to the large change in volume, fragmentation and disordering 
between the two phases can develop. In combination with a rigid lattice, this lead to a large 
hysteresis of phase transitions. This means that phases created under high pressure high- 
temperature conditions can persist inside the sample under standard condition ||T5[. Two 
other phases, the rhombohedral(r-BN) and the wurtzite (w-BN) are stable only at very 
high pressure (P~10GPa) and will not be considered here. Structural and excitation 
properties of wurtzite structure has been studied elsewhere [^]. At ambient pressures an 
interplay between h-BN and c-BN domains only should be expected Since 1979 

c-BN chemical vapor deposited ffims have been realized |jl9l, but the production of pure 
c-BN thin ffims (either by chemical or physical vapor deposition) remains a difficult task 



due to the formation, during the growth process, of unwilled h-BN domains [|18 . 

Extensive theoretical studies have been performed on the ground-state properties of 
BN. All calculations founded on density functional theory (DFT) and on the local densisy 
approximation (LDA) for the exchange-correlation potential agree in predicting that the 



c-BN structure has a lower energy than h-BN by about 0.06eV per atom pO| , |2T|p^ , |23[] . This 
result, which is con&med by our calculation, disagrees with an older calculation based on 
the orthogonal linear combination of atomic orbitals (OLCAO), where h-BN was found to 



be more stable than c-BN by 0.35 eV per atom [^. This is probably explained by a lack of 



convergence in the older calculations. Also the band structure properties of BN have been the 
subject of extensive thoretical work, performed within different schemes and approximations 
|l|gj|,|,|ig,gg,g3,0,0,g,0,0,0,g|,^,|g,0,0,||,g,||,g,g. On the experimental point 



of view many methods have been employed to explore the electronic structure. Among them 
there are methods like soft-x-ray emission spectroscopy (SXRES) P5| , |5P| and photoelectron 



spectroscopy(PES) [Q, which are both sensitive to occupied states, as well as near-edge 
x-ray absorption fine structure spectroscopy(NEXAFS) and energy loss near edge 

spectroscopy, which can be used to probe the empty bands(ELNES) In addition, 

electron energy loss spectroscopy performed in the plasmon region has proved to provide 
valuable information g|j4|g^ . 

The experimental optical functions of BN compounds are not as well known as those of 
the other group III nitrides. This is due in large part to the lack of high quality, single crystal 
samples. Consequently, most optical studies on BN have been performed on polycrystalline 
samples, with various amount of impurities. As one could expect, the results of optical 
studies performed on these materials vary greatly However a general consensus 



from the experimental data is that the minimum band gap is direct in the case of h-BN, 
and indirect in the case of c-BN. For h-BN a direct band gap of 5.2 ± 0.2eV associated 
to the transition Il3t,-Il2c has been extimated E3], while a value of 6.4 ± 0.5eV for the 



indirect minimum band gap in c-BN has been determined and associated to the Fis^- 
Xic transition [p5| , p3 . 



In the present paper we present ab-initio linear optical functions of c- and h-BN studied 
within DFT-LDA theory. No external parameter has been used to fit the experimental 
curves. Due to the importance of the two BN phases, and to the fact that in several 
experimental configurations they appear together in the same sample, a parallel study of both 
phases has been accomplished throughout the present paper. Complex dielectric functions, 
refractive index, refiectance spectra, optical conductivity, and energy loss functons for c-BN 
and h-BN have been calculated. We report our results in the energy range of interest to 
compare them with the existing theoretical and experimental literature. Calculated DFT- 



LDA energy-loss functions agree with experiments and with previous theoretical results, both 
for h-BN and c-BN. On the other hand, discrepancies arise between theoretical results and 
experiments in the imaginary part of the dielectric function for c-BN. Possible explanations of 
this issue are proposed and evaluated in detail. The imaginary part of the dielectric function 
in the h-BN case, instead, shows a reasonable agreement with both other theoretical results 
and experiments. 

This paper is organized as follows. In Section II we give the computational details of 
the calculation of the ground state properties for both phases. In Section III we describe 
the theoretical scheme used to evaluate the optical properties, and the part related to the 
calculation of excitation properties. In Sections IV and V we present results for the optical 
properties of c-BN and h-BN respectively. Finally, in Section VI conclusions are drawn, and 
the perspectives and issues opened by our work have been presented. 



II. COMPUTATIONAL DETAILS AND GROUND STATE CALCULATIONS 

Density functional calculations have been carried out within the local density approxima- 
tion (LDA) for the exchange and correlation functional [Q, using the Perdew and Zunger 



parametrization of the Ceperley and Alder results Kohn-Sham orbitals are ex- 

panded in a plane-wave basis set, with an energy cutoff of 55 Ry. 

Care has been used in constructing the ionic pseudopotentials, in order to avoid the 
occurrence of ghost states and to assure an optimal transferability Angular com- 

ponents up to / = 2 have been included. Separable, norm-conserving soft pseudopotentials 



have been generated within the scheme of Martins and Troullier ||5^ , with the following core 
radii (bohr): 1.59 (B, Sc?), 1.49 (N, 2p). For Boron, non linear core corrections (NLCC) to 
be used in the solid calculation, have been taken into account in generating pseudopotentials 



and pseudowavefunctions |56|. 



With the selected plane-waves cutoff, both the total energy and Kohn-Sham eigenvalues 
are converged to better than 0.1 eV for both c-BN and h-BN. For the cubic phase, the use of 
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10 special Chadi and Choen k points for charge integration in the first Brillouin Zone (BZ) 
||57| , is found to be sufficient to achieve a good accuracy for the computed total energy; for 
instance, total energy changes by less than 10 meV, and the fundamental band gap by less 
than 0.1 eV, when passing from 6 to 10 special points for BZ sampling. For the layered 
haxagonal case the use of 12 k points are sufficient for a full convergency of the total energy 
and eigenvalues. 

Hexagonal BN, turns out to be less stable than the cubic BN by about 0.050 eV/atom. 
Preliminary calculations in which "hard" norm-conserving pseudopotentials of the Bachelet, 



Hamann, Schliiter (BHS) form were used with energy cutoff of 150 Rydberg p5[ gave an 
energy difference of 0.078 eV/atom. Our results hence confirm the previous finding that 
c-BN is more stable than h-BN. No contributions to the total energy coming from zero point 



vibration of the lattices have been cosidered 18,20 



In Table(I) the structural parameters for the two phases are reported, in comparison 
with previous theoretical calculations and experiments. For c-BN our equilibrium lattice 
parameter underestimates the experimental one by few percent as usual for DFT-LDA, 
in accordance with previous theoretical results. The bulk modulus falls very near to the 
experimental energy range P0J2^ . 



An underestimation of the same order is found for the two lattice parameters a and ^ of 
h-BN, as reported by other authors too p0| , |24 |. 



No large differences in equilibrium structural parameters arise when the harder BHS 
pseudopotentials and the larger cutoff are used [p5|] . 

III. OPTICAL PROPERTIES AND ELECTRONIC EXCITATIONS 

The experimental scenario about the electronic properties of Boron Nitride is character- 
ized by many measurements. For h-BN, however, several old optical and energy loss data 
have been very conflicting each other mostly because of the use of poor samples. Band 



structures calculated in the past were not accurate enough to be of real help EM. Recently, 
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careful measurements of inelastic-electron-scattering spectra on well characterized h-BN be- 
came available |^^, and accurate measurements of the linear optical properties of c-BN 
have been reported [^. Ellipsometry measurements have also been done to characterize 
noncrystalline thin films of BN |]^ , |58| , and optical reflectance spectra of c-BN have been 
recently published . 



We determine the optical properties of h-BN and c-BN by calculating the momentum 
matrix elements (MME) associated with dipole transitions at a large number of k-points 
in the BZ. The frequency-dependent imaginary part of the dielectric function is given by 

US 



^ \ <v,k\ \c,k> fS{Ec{k) - E^{k) ~huj) 

k 



(1) 



where v and c label the valence and conduction states associated with the energies E^^k), 
Ec{k) , V is the crystal volume, and <> is the matrix element of the momentum operator. 
Nonlocality effects in the ionic pseudopotentials are neglected in the evaluation of e2 
Eigenvalues and eigenfunctions appearing in Eq.(l), are those determined within the DFT- 
LDA scheme. As it appears from Eq.(l), local field effects are also neglected in the dielectric 



function calculation |5y]. Reference formulas for the macroscopic optical functions, namely 

m 

N{uj) =n{uj)+ik{uj) (2) 

for the complex refraction index, will be frequently used in the following. Its relation with 
the complex macroscopic dielectric function reads: 

N^{uj) = e{uj) (3) 

with 

eiiio) =n\uj)-e{uj) (4) 



and 

e2(cj) = 2n{uj)k{uj) . (5) 
The absorption coefficient reads (c is the speed of hght) 

ri{uj) = 2ujk{uj)/c (6) 
which can be directly connected with the imaginary part of the dielectric function 

r]{uj) = uje2{uj) /n{uj)c (7) 
and to the optical conductivity a 

a{uj) = ri{uj)n{uj)c/A'K (8) 

In the above formulas once t2 (or k{uj)) is calculated or measured, the correspondent real 
part ei (or n{ijj)) can be determined via the Kramers-Kroenig transformation [|6l|] . 

To give an overview of the optical properties of the two crystals, and in particular of 
the differences among the hexagonal and the cubic phase, the calculated interband optical 
conductivity a is reported in in Fig. |I| and in Fig. |^. In the case of h-BN, an average over the 
three crystallographies directions is taken. The Joint Density of States (JDOS) for the same 
samples are also shown. Both the absorption threshold and the oscillator strengths in h-BN 
are much smaller than in c-BN. In the two cases, our data fairly reproduce the ortogonalized 
LCAO ones by Xu and Ching [Q. In both crystals there is a very little resemblance 



between the optical conductivity and the JDOS curves, especially in the high frequency 
region, indicating a strong k dependence of the MME. Both phases show a fundamental 
indirect band gap, but in h-BN the indirect gap is smaller than in c-BN This 
leads to the lower absorption threshold in h-BN. Many structures can be directly related 
to measured ones, like the peak around 6 eV for h-BN which can be assigned to transition 
between bands in the in-plane directions. This correspondence reflects the strong anisotropy 
in the optical absorption of this material, which shows up even when macroscopic quantities 
are taken into account (Eq.(2)-Eq.(8)). 
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As is well known DFT-LDA band structures for semiconductors and insulators can- 
not reproduce the real (experimental) ones. DFT-LDA eigenvalues, when interpreted as 



quasiparticle (QP) energies show the so called band gap problem ||62|. The QP energies cor- 
respond only qualitatively to the DFT-LDA ones, mainly because the band gaps between 
conduction and valence bands show a systematic underestimation with respect to the ex- 
periments |l63i|64| , |65i . This problem can be solved, for semiconductors and insulators, within 
a self-energy scheme, called GW Within this method, the self-energy operator reads 
S = GW, where G is the one-electron Green function and W the screened Coulomb in- 
teraction in the system, fully taking into account the screening properties of the material 
||6^ , |66|j67| . The corresponding Dyson equation can be solved within first-order perturbation 
theory with respect to (S — V^c), where V^c is the DFT-LDA exchange and correlation poten- 
tial. Assuming the DFT-LDA eigenvalues and eigenfuntions as the zeroth order eigenvalues 
and eigenfuntions of the system respectively, the QP corrections for the Bloch state ipnk 
state read 11621: 



where S'^^^ and V^^"^ are the expectation values of the Coulomb Hole (COH), Screened 

nk ' nk xc,nk ^ \ / j 

Exchange (SEX) contributions to the self energy, and of the DFT-LDA exchange-correlation 
(XC) potential, respectively. [^] Ab-initio GW schemes are based on the full calculation of 
the screening function of the system (i.e. fully including local field and dynamical effects) 
starting from DFT-LDA eigenvalues and eigenfunctions. These kind of calculations for 
semiconductors and insulators have often led to a very good agreement between theory and 
experiment []62| , |67 . 



Self-energy corrections to the DFT-LDA band structure for semiconductors and insula- 
tors can also be performed, with very good results, also using a model dielectric function to 
mimic the real screening of the system [^,0,^ . One can start for a perturbative formula 



of the form: [3S 
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nk 



i + di:^^,iE)/dE\ 

where Ti'^^^{E^^}) is the zeroth order term of the energy expansion of the dynamical part 

Tin 'flifC 

of S, and only the static parts of the SEX and COH self-energy are present p^ , |68[] . Cor- 
rections to the DFT-LDA spectra calculated from Eq.(lO) with the use of model dielectric 
functions enabled both us and other authors to obtain very positive results, when compared 
with experiments and ab-initio GW ones, in several semiconductors and insulators |]69| , |84| . 
Moreover, the use of model screening functions enables one to perform GW calculation in 
a more efficient way, reducing the computational effort of at least one order of magnitude 
in CPU time [^,0]. Details concerning the evaluation of the different coefficient appearing 



in Eq.(lO) and how do they compare with the ab-initio calculated ones are given elsewhere 

MM- 

Either ab-initio or simplified GW schemes have been employed to calculate the QP en- 
ergies of h-BN and c-BN [^,^,^,[7^]. All theoretical results improve the agreement with 
experimental data. Tipically, LDA transition energies suffer of an underset imat ion of about 



3eV for the cubic phase, and 2eV for the hexagonal one ||2^. In Tab. (II) and Tab. (Ill) for 
c-BN and h-BN respectively we report the results of previous calculations and experiments, 
namely the fundamental gaps and the main transition energies, in comparison with present 
results. In c-BN the present transition energies have been calculated at the corresponding 
theoretical lattice costant while in ref. [^, (as explicitely stated by the authors) the ex- 
perimental lattice parameter has been used. Moreover, the use of an efficient method in a 
large gap system induces a slight overestimation of the above energies [0,^ . As it is clear 



from the data, minor differences arise between the ab-initio GW results and the efficient 
GW ones. 

In Tab. (Ill), the value marked with "a" correspond to a transition between a point Ti 
(along the T — K direction near K, and near in energy to H^,) and Mc. Similarly the term 
with "b" is the transition energy between K^^-K^c, which is near in energy to the transition 

H3t,-H2c p3 |. 
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IV. CUBIC BN 



In Fig. ^we report the calculated dielectric function for c-BN obtained within DFT-LDA. 
The static dielectric constant is found to be of 5.45, larger than the experimental one (4.45) 
||7l| . This fact must be ascribed to absence of local fields effects in the calculation of €2 |l36| , |60| . 
Experimental curves obtained from recent optical reflectance measurements in the energy 
range 5 -i- 25 eV with syncrotron radiation by Osaka e coworkers [^] are also reported. Two 
different samples have been used by them, namely a sintered c-BN plate and c-BN thin film. 
The c-BN sintered plate produced at high pressure and temperature, offered by Sumimoto 
Electric Industry (Japan), had dimension of 5x5x0. 5mm^ and was characterized as cubic by 
using X-ray diffraction and Raman Scattering [^]. The film, on the other hand, was grown 
using chemical vapor deposition (CVD) on a silicon substrate, with negative self-bias. c-BN 
phase, syntesized from B2HQ {At diluition) and N2 gas mixtures, had a the grain size is 
100-200 A, controlled by Transmission Electron Microscopy (TEM) images. 

In Fig. the calculated reflectivity of c-BN is displayed together with the experimental 



data taken from ref. Both the spectra show peaks at 11.7 and 14.0eV and a broad 

structure around 18eV. The authors of ref. |^6|| claim that the first two peaks correspond 



to the El and E2 peaks of the zinc-blende- type semiconductors [£T|. Our theoretical curve 
shows two structures centered at 12.8 eV and 14.4 eV but no evidence of higher energy peaks 
results. To enter in further detail with other theoretical data and experiments, we report in 
Fig. I the imaginary part of 62, in comparison with two experimental curves by Osaka et al. 
46| One of these curves has been obtained from a crystalline sample, and the other from a 



film grown by CVD technique. In Fig. ^ we also compare our curve with that calculated by 
Xu and Ching ||23], and with a LMTO one calculated at our theoretical lattice parameter 



72|. It is significant to discuss these figures in detail. The experimental data show major 



structures A and B at 9.05 and 11.7 eV, and shoulders at 13.2eV, 16.7eV (C and D) for the 
crystalline sample. In the €2 in the c-BN film, the peak A and the shoulders C and D are not 
present. The spectrum of the c-BN thin film deposited by plasma CVD is similar to that 
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of the c-BN synthesized under high pressure and temperature. The authors claim that the 
differences within their spectra is a demonstration of the fact that the microcrystalhne c-BN 
does not grow epitaxially on sihcon substrate, and presents a rough surface. They do not 
mention the value of the absorption onset, but a value of 6.8 eV can be deduced from their 
figure of €2- To consider the correspondence with the electronic transitions, peak A has been 
associated to the transition Fis^, — Fisc, the transition X^^-X^c and X5t,-Xic, peak C to the 
transition L^y-L^c, and peak D at 16.7 eV to Lit,-Lic [|^. The present DFT-LDA reflectance 
spectrum reproduced qualitatively well the experimental curves, but the first peak is found 
at 12.8eV and the following at 15.2 eV. Fig. ^ also demonstrates that our results do 
agree with other theoretical calculations (moreover the theoretical curves reported show no 
remarkable difference in the range 25 — 40 eV which has not been showed in this figure). 
However, the absorption onset obtained within DFT-LDA is at 8.9 eV, i.e. much higher 
that the experimental one. This issue turns out to be rather singular in sight to the the 
fact that usually , due to the above mentioned band gap problem, the DFT-LDA calculated 
absorption onset should fall at /ower energies than the experimental ones P5| , |5B| , P^ , |52| . In 



our case the LDA theoretical curve is blueshifted by more than 2eV with respect to the 
experiment. A similar overstimate is found also by other theoretical calculations P^ , |5B| , |55| . 
The interpretation of Christensen and Gorczyca |^ for 62 is different (away from treshold) 



from that proposed by Osaka e coworkers , Xu and Ching 0] and from our interpretation. 



Infact the authors of ref. |^ associate the A peak to the transition 3 to 5 at the 1/2TU 
point, the major peak B to the 4-5 transition along S and L point, the C peak to the 4-6 
transitions at L, and the D peak to transition 4-6 at A and at S'. The onset is ascribed to 
direct transition at F. Xu and Ching found A, B and D structures respectively at 10.7eV, 
12.6eV and 15.6 eV, but no evidence for the C shoulder in their curve results Also 



their data show a rigid blueshift with respect to the experimental curve by about 1.7eV. 
Our results rensemble those of Xu and Ching. Please notice that their name assignements 
for the spectrum does not match the present ones. Our onset is found at 9.05eV, while the 
structures A, B, C and D are found respectively at 9.7 eV, 12.5 eV, 14.9 eV, with small 
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evidence of peak D. Our assignement to band transitions within DFT-LDA agrees with that 
proposed in Ref. which has also been confirmed by Tsay and coworkers ||27|| . 



As far as we know the only other existing experimental data for the minimal direct gap 
of c-BN are those by Philipps and Taft [Q, which yield a value of 14.5 eV. They studied 
diamond in the vacuum ultraviolet, and they reported data about a c-BN sample displaying 
structures around 9 and 10 eV, and a peak near 14.5 eV The BN data were reported as 
marginal ones by the authors, to support the evidence for a larger direct band gap in c-BN 
than in diamond. Their value of 14.5eV has been considered as the experimental direct 



band gap of this material by other authors in the literature [p3| , [71| . On the other hand, 
these data could be interpreted in good agreement with our DFT-LDA spectra. Infact, 
one could consider the first structure as the display of the onset (at 9eV in our DFT-LDA 
calculation) and the major peak as the maximum (at 12.5 eV in our DFT-LDA calculation) 
(see Fig.s(2)-(4)). 

Let us now address in detail the main mismatch between theoretical results (including 
ours) and the experiments, concerning the onset. One of the points to be considerd is 
the problem of the correct value of the lattice parameter to be used in optical properties 
calculations. This issue has been raised for BN by experimental and theoretical works 



ll,!!. In Fig.(0) we reported different DFT-LDA curves for 62 calculated at different 



lattice parameters around the equilibrium one. It is evident that the onset of the theoretical 
absorption meets the experimental one only when the lattice constant is expandend, with 
respect to the equilibrium one, by 5%. This value appear to be too large, both with respect to 
the experimental one, and in sight of the fact that: I) usually the DFT-LDA understimation 
of the lattice costant is about (1%) [^,0, and II) usually the DFT-LDA understimates the 



absorption onset. 

Another possibility is that the sample under experimental study has not been carefully 
characterized, namely concerning different possibilities of disturbances, as thermal ones, or 
the presence of (oxide) overlayers PB| , ^ , which we will not address here but which cannot 



be excluded in principle. In the following, we consider instead the possibility of h-BN 
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contamination within the c-BN sample. This is possible due to the fact that in normal 
conditions h-BN domains may be found in c-BN samples 

In Fig. ^ and Fig. |^ we also report the real and imaginary part of the refractive index, 
A^(ti;) and the reflectance spectrum in comparison with the experimental data of Miyata and 
coworkers The experimental data have been obtained from reflectance measurements 
over a photon energy range between 2 23eV, and from transmittance data in the range 
2 7eV. A single crystal of c-BN (5mm^ area x 0.16mm thick) has been used. Both the 
onset of the absorption coefficient, or the imaginary part of the refractive index, indicate a 
gap of 6.1 ±0.5el^. Our results are only in qualitative agreement with these measeurements. 

We can argue that part of the error comes from the normalization imposed to the ex- 
perimental curves. Infact, the absolute value of reflectance and transmittance have been 
obtained by the authors of ref. |^5| by imposing that refractive index aX u = 2.10 eV equal 
to 2.117, the value determined by an indipendent experiment by Gielisser et al. Within 
the present work instead, no adjustable parameter has been used to fit the experimental 
data. 

In Fig. [ly we report the calculated DFT-LDA energy loss function for cubic BN , namely 
/m[— l/e(g, u)], in comparison with the experiments. We compare our results with Electron 
Energy Loss (EEL) mesaurements of McKenzie et al. [^. The dominant maximum in the 
EEL spectrum is due to the plasmon excitation, a longitudinal oscillation of the valence 
electrons as a whole against the cores, classically occurring at frequency = Ne^/m were 
N represents the density of the valence electrons of the sample. Due to the different density 
of h-BN and c-BN, EEL spectra have been used to discriminate between the two phases 
within the same specimen The measured Up for h-BN is 25.5eV and for the c-BN is 
28.5eV. An absorption threshold for the cubic phase at 9eV has been deduced by the authors 
of ref. Our DFT-LDA calculation yields a first peak at 30eV, which is in reasonable 

agreement with the experimental one (28.5 eV) and with that determined by Xu and Ching 
(28eV) 
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V. LAYERED HEXAGONAL BN 



For this material conflicting experimental and theoretical values of the fundamental band 
gap can be found in the literature. Tarrio e Schatterly presented EELS spectra between 
-r-60 eV, with a plasmon peak at 26.4eV. From their data, a direct gap of 5.9 ±0.2 eV was 
inferred ^7^. Catellani and coworkers using a linear augmented plan wave computational 
scheme (FLAPW-LDA), claimed that h-BN has a minimal indirect gap at the H^y — Mic 
transition, corresponding to 3.9 eV, and a minimal direct gap at 4.3 eV, corresponding to 
the H^^ — transition . Hoffman and coworkers, by optical reflectivity measurements 
in the range 0.045 -r- 10 eV, found a direct gap of 5.2 ± 0.2 eV, deducing that it correspond 
to direct transitions at H Park, Terakura, Hamada, using a FLAPW (LDA) scheme, 
found a minimal direct gap, and gave for the — transition energy a value of 4.5 eV 
||30| . Last but not least, the work of Suhr and coworkers and |3^ of Cappellini et al. found, 
within DFT-GW, that the material has an indirect minimal gap (see section III). For h-BN 
we start the discussion from the EEL function reported in FigjTl|. In this figure, we compare 
the calculated DFT-LDA curve, averaged over the three crystallographies axis, with the 



curve calculated by Xu and Ching and with the experimental one by McKenzie |[43|| . 
We recognize three major structures. In the theoretical spectrum of ref. [Q, a first peak(A) 
around 7eV can be found, a second peak (B) is at 12eV , and a major structure (C) appears 
at 24 eV . The latter is the bulk plasmon peak, to be compared with the experimental value 



of 25.5 eV In Fig. ^ we report the calculated imaginary part of the dielectric function 
for x-y and z components. Our results compare well also in this case with the LCAO ones 
by Xu and Ching. In Fig. O we plot the imaginary part of the dielectric function averaged 



over the three crystallographies axis. In particular after Fig. |I2| and Fig. |I3|, we agree with 
the conclusions of Xu and Ching that peak B in the EEL function of FigjTT|, comes from 
the component of €2 parallel to the c axis of the hexagonal crystal, while peak A comes 
form the perpendicular components of €2- These statements may be easily confirmed by 
considering experimental curves at different scattering momentum transfer (in plane versus 
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c-axis data) [^,0. The major peak C contains contributions from both the €2 components 
perpendicular and parallel to the c-axis. In Fig. |1^ we report the imaginary and real part of 
the calculated refractive index. Again, a qualitative difference with respect with the cubic 
case is observable. Infact, although the imaginary part of the dielectric function has been 
averaged over the three crystallographies directions, strong differences arises with respect to 
the cubic case (see Fig. (8)). 

The knowledge of the optical properties of h-BN can help us to understand the mismatch 
between theory and experiments for the €2 of c-BN. Infact, as it was recently shown |]74 



c-BN films display an hexagonal-like top layer, and in presence of disturbances the amount 
of disorder in both types of films increases significantly, leading to the transformation of 
the cubic phase to the hexagonal-like material. Moreover, it has been demonstrated that 
due to the large hysteresis of the hexagonal-cubic phase transitions, h-BN domains may 
continue to exist in c-BN sample . It can hence be argued that the experimental samples 



of c-BN are likely to contain impurities due to h-BN domains. Assuming that the sample 
used in the experimental work by Osaka and coworkers contained h-BN domains, we 



can work within the effective crystal approximation, i.e. assuming that the domains are 
homogeneous and isotropic: in this case an imaginary dielectric function corresponding to 
a linear combination of the two pure forms can be expected for the merged system [f75[] : 



e^H = xe^^-^^'icu) + (1 - x)6r^^(^) (11) 



with X going from to 0.5. In Fig. ^ we report the results obtained within this scheme 
always at the DFT-LDA level. By increasing the concentration of h-phase, a peak at 5 eV due 
to h-BN presence does grow, while the peak (for which the cubic phase is responsable) around 
12eV, due to the cubic phase decreases. Even if the main peak comes in better agreement 
with the experiments in intensity, the peak around 5eV remains significantly red shifted 
with respect with the experimental onset. To address more deeply this issue, we calculate 
self-energy corrections to the spectra of both phases within the GW approximation . 



The results are reported in Fig. [T^. The GW corrected spectrum matches the experimental 
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treshold better than the DFT-LDA one, but the major peak now misses completely the main 
experimental structure at 11.7eV. 

Another possible explanation to the mismatch between theory an experiment for the 
cubic phase is the contribution of higher-order effects in the spectra, like those due to the 



electron-hole interaction which we have neglected so far. If we consider, in Fig. ITBI, the case 
with x=0 (pure c-BN with GW self-energy corrections), the GW theoretical onset falls at 
11.6eV (see section 1), more than 5eV higher than the experimental one (~ 5.8eV by Osaka). 
Hence the inclusion of self-energy effects in the DFT-LDA spectrum worsens the comparison 
with experiments in the present case. An effect of this type is typical for systems where the 
excitonic effects are large. The excitonic binding energies infact, can range from few meV 
(like in bulk GaAs and Si) to much larger values (of the order of the eV) in insulating oxides 
||76| , [77|J78| . The effects of self-energy and excitonic corrections to the DFT-LDA €2 curve 
roughly correspond to shifts in opposite directions |176|J78| . Referring to the onset of 62, one 
might obtain a rough estimate of the excitonic binding energy from the mismatch between 
GW corrected onset and the experimental one. In the c-BN value of more than 

5eV would occur. This value seems to be too high, even in comparison with the strongest 
exciton binding energies found in systems like oxides, where, due to the weak screening, 



the electron and the hole can be bound by an energy binding of the order of leV ||76|. A 
theoretical evaluation based on Wannier functions for c-BN leads to a value of the this order 



of magnitude [Q. To fully address this problem from the theoretical point of view, one 
should solve the Bethe-Salpeter equation for the optical response function including two 



particle effects [|76| , [77| , |78| . On the other hand, one should look for more refined experimental 
mesaurements of the linear optical properties of c-BN. In summary c-BN and h-BN mixing 
within the sample, self-energy and excitonic effects, the presence of impurites, the lack fo 
pure crystalline samples, and lattice parameter mismatch, can in principle all be sources of 
the disagreement between theoretical results and the available spectra for c-BN. 
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VI. CONCLUSIONS 



Linear optical functions of cubic and hexagonal BN have been studied within DFT-LDA. 
Calculated energy-loss functions compare well with experiments and with previous theoret- 
ical results, both for h-BN and c-BN. Discrepancies arise between theoretical results and 
experiments in the imaginary part of the dielectric function for c-BN. Possible explanations 
to this issue are proposed and evaluated; lattice constant variations, h-BN contamination in 
c-BN samples, and self-energy effects. On the other hand our DFT-LDA results show a rea- 
sonable agreement with other theoretical outcomes and with experiments for the imaginary 
part of the dielectric function in the case of h-BN. More refined measurements and calcula- 
tions are needed to fully address the mismatch between theory and experiment for the cubic 
case. It is in our programs to go further in direction of a refinement of the calculations. 
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TABLES 



c-BN 


Workl 


Work2 


Theol 


Theo2 


Exp. 


ao{a.u.) 


6.754 


6.771 


6.759 


6.833 


6.833 


Bo{Mhar) 


4.01 


3.52 


3.97 


3.70 


3.69-4.65 


h-BN 












ao{a.u.) 


4.698 


4.68 


4.50 


4.71 


4.72 


c/a 


2.608 


2.6068 


2.608 


2.670 


2.664 


Bo{Mhar) 


2.68 


2.65 


2.61 


3.35 





TABLE L Calculated structural properties of zincblende and hexagonal BN (lattice constant 
and bulk modulus) compared with other theoretical and experimental results. First column: 
present calculations; second column: theoretical results by Cappellini et. al. |25|, Theol: re- 



sults by Furthmiiller et al. [20|, Theo2: results by Xu e Ching p3, and in the last column the 



experimental values from Ref. |25| 



c-BN 


GWl 


GWlbis 


GWf 


Exp 


Tis'i, — Xic 


7.28 


6.95 


6.3 


6.4 ±0.5 


Lis't, — Tie 


11.79 


11.46 


11.4 


14.5 



TABLE IL Excitation energies for c-BN. In the first column, the present theoretical values here 



calculated (GWl), GWlbis is taken from Ref. GWf from Ref. The experimental results 
for the indirect gap are from Ref. [^^, and those for the direct one are from Ref. j^. 
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h-BN 


GWl 


GWlbis 


GWf 


Exp 




6.39 


6.04 


5.4" 


5.2 ±0.2 


H3V — H2c 


6.76 


6.66 


6.33^ 





TABLE III. Excitation energies for h-BN . The first column contains the values calculated here 



(GWl), GWlbis is taken from Ref. [25|, GWf from Ref. |10| and the experimental values from Ref. 
For the meaning of the a and b marks, see text. 
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FIGURES 
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FIG. 1. Calculated optical conductivity a of c-BN (solid line) vs. the theoretical result of Ref. 
|2| (dashed line). 
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FIG. 2. Optical conductivity function of h-BN (solid line) vs. previous theoretical result 



(dashed hue) after Ref. ||2J. 
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FIG. 3. Computed real and imaginary parts of the dielectric function for c-BN. ei is the solid 
line, 62 the dashed one. 
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FIG. 4. Calculated reflectivity (solid line) of c-BN vs. experimental results in arbitrary units 
(dashed lines), taken from Ref. |Q. 
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FIG. 5. Calculated imaginary part of dielectric function of c-BN (solid line) vs. the experimen- 
tal results of Ref. ||4^ ( dashed lines). 
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FIG. 6. Calculated imaginary part of dielectric function of c-BN (solid line) vs. LMTO result 
after Ref. [f| (dotted line) and LCAO result (dashed line) of Ref. ||. 



32 



24.0 



18.0 



^2 12.0 




12 16 20 

ENERGY (eV) 



FIG. 7. Calculated imaginary part of dielectric function of c-BN at experimental (dotted line) 
, theoretical (bold solid line) and theoretical plus 5% (solid line) lattice constant. Experimental 



results (dashed lines) after Ref. |^ are plotted for reference. 
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FIG. 8. Calculated real part of refracive index of c-BN (solid line) vs. the experimental results 
(dot-dashed line after Ref. ||4^) , and previous theoretical result (dashed line after Ref. [p4| ). 
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FIG. 9. Calculated imaginary part of refracive index of c-BN (solid line) vs. experimental 



results (dashed line) after Ref. p5 |. 



35 



0.30 



g 

O 

z 

3 
u. 

<J) 

(0 

o 

>- 
o 

DC 
lU 

z 

HI 



0.20 



0.10 



0.00 




10.0 20.0 30.0 

ENERGY (eV) 



40.0 



FIG. 10. Energy-loss function of c-BN (solid line) vs. the experimental result (dot-dashed hne) 



after Ref. [43| and previous theoretical result (dashed line) of Ref. |24|. 
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FIG. 11. Energy-loss function of h-BN (solid line) vs. the experimental result (dot-dashed hne) 



of Ref. 1 43] and previous theoretical result (dashed line) after Ref. 
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FIG. 12. First panel: imaginary part of the dielectric function of h-BN (solid line) averaged on 
X and y directions vs. the same quantity after Ref. ||2^ (dashed hne). Second panel: imaginary 
part of the dielectric function of h-BN (sohd hne) on z direction vs. the same quantity after Ref. 
H (dashed hne). 
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FIG. 13. Imaginary part of dielectric function of h-BN (solid line) vs. previous theoretical 
result of Ref. |^] (dashed line) averaged over the three crystallographic directions. 
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FIG. 14. Calculated refractive index of h-BN: real (solid line) and imaginary part (dotted line). 
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FIG. 15. Imaginary part of dielectric function for merged c-BN and h-BN (solid lines). The 
experimental results of Ref. are also given (dashed lines). Pure cubic €2 is showed with bold 
solid line. 
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FIG. 16. Imaginary part of dielectric function with GW corrections (see text) for merged c-BN 
and h-BN (solid lines). Experimental results of Ref. |^] are also given (dashed lines). Pure cubic 
€2 is showed with bold solid line. 
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